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determination of the reliability of the net with respect to both 
complete communication and k-terminal communication. The effect on this 
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1 . 



Introduction. 



Any system that provides for communication, be it the transmission 
of information or the transportation of commodities, is subject to pos- 
sible failure. The loss of security, the malfunction of a radio trans- 
mitter, and the destruction of a bridge are examples of this failure. 
Linear graph theory, as described in texts such as [5], has been ap- 
plied to measure the reliability of the system given some measure of the 
reliability of the parts. Fu and Yau, [4], have indicated a value for 
the probability that two stations can communicate; and a digital computer 
algorithm for this has been given, [6]. The probability that all pairs 
of stations can communicate within a given time interval has been studied 
by Chan, [2], and Fu, [3]. 

In this work the application of the k-terminal reliability function, 
to be defined in Section 3, is presented, and a digital computer pro- 
gram is included. In addition, the application of the power of this 
function without computer assistance is demonstrated. A formula is given 
to find the value of the 2-terminal reliability function without enumer- 
ating the Table of Combinations. 
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2 . 



Statement of the Problem. 



A communication net may be represented by a linear graph in which 
the edges of the graph represent the communication links, and the nodes 
of the graph represent the stations of the communication net. 

A weight can be assigned to each edge of the graph to indicate the 
probability that the particular communication link represented by the 
edge will operate successfully. This weight is called the reliability 
of that edge. This probability is not a conditional probability. Each 
edge is assumed to be associated with an independent random variable, 
having only two possible states: normal operation, and not operational. 

The communication net can be analyzed by studying the topology of 
the graph. Topological methods can be applied on several different 
problems. For example, if one link of a communication net has been lost 
(destroyed), it must be decided whether or not to replace this link with 
one of the other existing links in the net. The information needed to 
make this decision must include the relative importance of each link, 
some measure of the quality of the branch (e.g., its reliability), and 
the role of the branch in the net. The analysis of the data must in- 
clude the degree of importance attached to the overall net reliability 
as contrasted to the reliability of communication between a pair or a 
group of terminals. 

Several investigators, [3], [4], [6], have proposed quantitative 
measures for the reliability of a communication net. Investigations by 
Chan [2], have led to a result which uses the reliability function as a 
k-terminal reliability function. This function may be applied to the 
reliability of communication between any k-number of nodes. The topo- 
logical methods which are used to generate the function are suitable 
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for use on a digital computer. The large storage requirement is the 
limiting factor in their application. 
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3. Properties of Reliability Functions. 

3.1 2-terminal Reliability Function. 

The reliability function is defined in reference [2). Some of 

the most important definitions are included here for ready reference and 

emphasis. Using the notation of [2], there is a value of reliability, 

p. , p„,...p , assigned to each of the edges of the graph, e , e ,...e , 

1 l e i 2 K 

...e . Each branch of the net is assumed to be associated with an in- 
e 

dependent random variable. If the branch, b^ , is operational, the edge, 

e^, of the graph assumes the value p^, which is the probability that 

edge, e,, will operate normally. If the branch is not operational, the 

edge, e^, assumes the value which is the probability that the edge, 

e , will not be operational. Although these reliabilities are probabili 
i 

ties, they are considered not to be conditional probabilities. That is, 
we have the following relations: 

P ± +J i = l (3.1) 

Pr (A and B) = Pr(A)Pr(B) (3.2) 



Definition 3.1 : 2-terminal Reliability Function : 

A reliability function expressed in terms of the edge reliabilities of 
a network is a 2-terminal reliability function if it represents the 

j| if 

reliability between two distinct terminals (i.e., the probability of 

communication between these two terminals) . 

Let q represent the 2-terminal reliability function between 
x y y 

terminals x and Then we have 



x >y 



" m lPl p 2 . 



+ m 2 P lP2 ... p e _ lPe 



+ m 2 e Pl p 2 . 



(3.3) 
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where m^, n^* • • . .m e are either 1 or 0 depending on whether the com- 
bination of states of the edges provide a communication path between ter- 
minals x and 2 * If the path exists, then m^, = 1; if not, m^ = 0, When 
the function is expressed as in equation (3.3), it is said to be in its 
canonical form . 

Definition 3,2 : Path Product . 

A path set is a set of branches of a communication path between a 
pair of terminals, A path product is the product of the variables 
associated with the elements in a path set. 

Definition 3,3 : Path Product Term or Primary Path Product Term , 

A canonical term is a path product term if 

(1) its unbarred variables are those of a path product; and 

(2) all other variables in the term are barred variables. 

Definition 3.4 : Secondary Path Product Term . 

A canonical term is a secondary path product term if its unbarred 
variables correspond to those elements forming a path set plus one or 
more not in the path set. 

Definition 3.5 : Table of Combinations . 

The collection of all nonzero canonical terms of a 2-terminal 

reliability function is called the Table of Combinations. That is, 

all the terms in which m. takes on the value 1. 

x 

Definition 3.6 : k-edged Path . 

A k-edged path is a path of k number of edges or branches which 
constitute a communication path between a pair of terminals. 



11 



Theorem 3.1: 



Let P. be an A-edged path between terminals x and £, in a graph 

A 

containing e number of edges. The Table of Combinations for the 2- 
terminal reliability function, q , will include a primary path product 

x >y 

f 6 —A ) 

term for path and 2 V -1 secondary path terms. 

Proof : 

From definition (3.3) and (3.5), the primary path product term 

(6*A) 

must be present. There are 2 V ' possible combinations of the (e-A) 

variables not specified by the primary term. However, definition (3.4) 

requires the number of unbarred variables to be at least (A+l). There- 

(g -A) 

fore, the maximum possible number of combinations is 2 -1. 

If a graph contains e number of edges and n number of vertices, then 

definition (3.1) indicates that there will be 2 number of terms in the 

2-terminal reliability function. Theorem (3.1) indicates that one path 

( e -A) 

of A number of edges will produce 2 terms in the Table of Combina- 

tions. This is the same as saying that, from definition (3.1), nu will 
( e -A) 

be nonzero in 2 terms of the 2-terminal reliability function. 

Lemma 3.1 : 

If no path exists between two terminals x and y, then the 2-ter- 
minal reliability function will have no nonzero terms. 

Proof : 

By definition (3.1) all the m^ values will be zero, since there is 
no path between terminals x and y. There will be no entries in the Table 
of Combinations. 
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Theorem 3.2: 



If a graph contains e number of edges and n number of vertices, and 

if an A-edged path is the only path between vertices x and y, then there 
(g -A) 

will be 2 number of terms in the 2-terminal reliability function, 

q 

x ,y 

Proof : 

( 6 * A) 

By theorem (3.1) there will be 2 number of terms for the A- 

edged path. Since there are no more paths, lemma (3.1) shows there will 
be no additional terms. Hence, the theorem. 

Theorem 3.3 : 

Let be an A-edged path and P^ a B-edged path, each between ter- 
minals x and £, in a graph containing _e number of edges. The primary 
path product term of path set A can not be a secondary path product term 
for path set B in q 

x,y 

Proof : 

The unbarred elements in the secondary path product terms must ap- 
pear as barred elements in the primary path product term of the path set 
A, by definition (3.3). Hence, the theorem. 

Theorem 3.4 : 

Let P A and P g be path sets between terminals x and y_. Then at 
least one secondary path product term for path set A will be identical 
to one for path set B. 

Proof : 

From definition (3.4), the collection of secondary path product terms 
for any path set will include the case where all variables are unbarred. 
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This will be common to all collections of secondary path product terms. 
Hence, there will always be this term present. 



Example 3.1 : 




Figure 3.1 



Path Set A 
Path Set B 




2-terminal reliability function, q .expressed in canonical form: 

x »y 



ra l p l p 2 p 3 


+ 




(7) (10) 




m 2 P l P 2 P 3 


+ 




(8) 




ra 3 P l P 2 p 3 


+ 




(9) 




m 4 P l P 2 P 3 




(2) 




(ID 


ra 5 P l P 2 P 3 


+ (1) 








m 6 P l P 2 P 3 


+ 


(3) 




(12) 


m 7 P l P 2 P 3 


4- 


(4) 




(13) 


m 8 P l P 2 P 3 




(5) 


(6) 


(14) 



(1) primary path product term for path set A 

(2) primary path product term for path set B 

(3) (4) (5) secondary path product terms for path set A 

(6) secondary path product term for path set B 
observe that (5) and (6) are identical 

(7) (8) (9) m^ terms are 0 because terms do not contain a path 

between terminals x and £ 

(10) this term will always be zero 

(11) (12) (13) a tree of the graph is made up of the unbarred 

terms 
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(14) a tree plus a chord of the graph is made up by the un- 
barred terms 

3.2 Multi-terminal Reliability Function. 

Reliability functions, which are expressed in terms of the edge 
reliabilities, may be derived in order to represent the reliability of 
transmission between more than two terminals. Thinking of a communica- 
tion net, one can envision a classification of the terminals into the 
categories of transmitter, receiver, or both, and repeaters. The re- 
peaters may be looked upon as intermediate steps in the communication 
paths. All of this may enter into an evaluation of the reliability of 
simultaneous communication between all pairs of terminals. 

Definition 3.7 : Complete Communication . 

The event that every pair of terminals x and £ in a net can communi 
cate with each other simultaneously is complete communication; the simul 
taneous communication between all pairs of terminals of the net. 

Definition 3.8 : N-terminal Reliability . 

The probability of complete communication in a net with N-terminals 

is called the M N-terminal reliability" of the net and is denoted by q . 

n 

Definition 3.9 : k-terminal Reliability . 

In a set of n number of terminals, consider k number of terminals 
of transmission, where k < n. Then, the probability of complete com- 
munication among the set of Jc number of terminals is called the k-ter- 
minal reliability of the net and is denoted by q . 

K 

From definitions (3.7) and (3.9), a net may have several differ- 
ent values of k-terminal reliability depending on which terminals are 
included in the set of k number of terminals of transmission. The 
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choice of terminals included in this set is of fundamental importance 
in the interpretation of the analysis or data collection. 




For a given net, q^ is found in the following way: first determine 

for each pair of vertices. There will be %(n)(n-l) for a net with 



n number of terminals. Then sum all the canonical terms which are common 

to all q functions. To find q , find q for each pair of terminals 
x,y k x,y 

in the set, K. There will be %(k)(k-l) number. Then sum all the canoni- 
cal terms common to all these q functions. 

x,y 

Another method is first to determine a complete set of trees for the 
net. Complete communication is possible if and only if the unbarred vari- 
ables of a canonical term include those corresponding to the branches of 
a tree. By taking the sum of all canonical terms that contain a tree, 
or a tree plus chords, made up from the unbarred variables, q^ is deter- 
mined. If the terms which are summed contain a subtree of the graph, or 
a subtree plus chords, containing all the vertices of the set K, the q^ 
is determined. This method can be applied on the canonical form result- 
ing after the first q calculation. This can reduce the labor required 

x,y 

to find q^ or q^ provided all the trees of the graph are known. 




Path set P c = 

Path sat P D - 

Path set P £ = $ e jj 



X * 

Figure 3.1 repeated 



P F " ^®l e 2 } 



2 



Path set 



m 5 P l P 2 P 3 + 



m 6 P l P 2 P 3 + ^ ^ m 6 P l P 2 P 3 + ^ (8 ^ 

m 7 P 1 P 2 P 3 + (4) (9) m 7 p l p 2 P 3 + ^ ^ 

m 8 P l P 2 P 3 (5 > (1 °^ m 8 P l P 2 P 3 ^ (10 ^ 



(1) Primary path product term for path set C and E. 

(2) Primary path product term for path set D and F. 

(3) (4) (5) Secondary path product term for path set C and E. 

(6) Secondary path product term for path set D and F. 

(7) (8) (9) (10) make up a tree, or a tree plus a chord, of 

the graph, made up of the unbarred terms. 

Example 3.3 




Figure 3.1 repeated 

The 3 trees of the graph are: 



Number of trees, t 




-1 

2 



= 3 



^®l e 2^ ’ { e 2®3^ ’ [ e i e 3^ 

Using the results of examples (3.2) and (3.1), we have 

% ■ Vi'* + 

m 6 p l p 2 p 3 + 
ra 7 P l p 2 P 3 + 
m 8 p l p 2 p 3 

by the process of summing the canonical terms in common. From example 
(3.1) alone, or (3.2) part (a) alone or part (b) alone, the same results 
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could have been obtained using the method of trees. 

In this simple example, for any set of terminals of transmission 
less than the given set, the problem reduces to one of the examples, 
since it becomes a 2-terminal reliability problem. 



3.2 Some Characteristics of the Path Product Terms. 

A graph contains 5 edges. A matrix can be constructed in which the 
elments are the number of path product terms that contain U number of un 
barred variables compared to the L number of edges in the path. 

Number of unbarred variables 



Path Length 
4 
3 
2 
1 



1 2 3 4 5 

0 0 0 1 1 

0 0 12 1 

0 13 3 1 

1 4 6 4 1 



The number of path product terms 
This array is recognized as the start of Pascal's Triangle. See 
reference [1]. The nonzero matrix elements are the so called binary 
coefficients, defined by the formula 




n! 



(3.4) 



j -'(n-j).' 

where, in this case, n = (e-A), and A is the path length, and j = 0, 
1,2,..., (e-A). 

If all the branches of the net have the same reliability, p, then 
the collection of all path product terms resulting from an A-edged path, 
in a graph with e edges, can be written as follows: 
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( e - A ) pV*- a) + (•;*) p c** i »-c^-i> + M 



e-A-1) + fe-A] p (A+2)-(e-A-2) + 



<“■) 



e-U-Cl) +(«-*]„« 






e-A 






i - 0 



(3.5) 



where K(i) is the collection of the binary coefficients defined by the 
formula 




n: 



i ! (n-i ) ! 

Examples (3.1) and (3.2) illustrate the generation of the primary and 
secondary path product terms. Theorem (3.4) points out that there will 
always be at least one of the secondary path product terms common to 
the expansions generated from two or more path sets is a 2-terminal re- 
liability function. Theorem (3.3) proves that the primary path product 
term is never a secondary path product term for another path in the 2- 
terminal reliability function; that is, it is not one of the duplicate 
terms. The exact number of unique terms added to this function is a 
problem in conditional probability. The following group of theorems is 
fundamental to the solution of this problem. 

The theorems and lemmae in this section are subjected to the follow- 
ing conditions: a graph contains e number of edges and n number of verti- 
ces. Path set, P^ is a 1-edged path set between terminals x and £, 

made up of branch Id, alone. 
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Theorem 3.5: 



No multi-edged path between terminals x and y may contain branch b. 
Proof : 

By definition (3.2) a path set is a set of branches of a communica- 
tion path between a pair of termunals. No additional branches are in- 
cluded in the path set. Hence, since branch b is a path by itself, any 
set that is made up of branch b and one or more additional branches can- 
not be a path set between terminals x and Definition (3.6) of a k- 
edged path allows the conclusion that no multi-edged path between ter- 
minals x and £ contain branch b. 

Lemma 3.2 : 

There can be no primary path product term in which branch b ap- 
pears as an unbarred variable in the 2-terminal reliability function, 

q , except path set P 
x , y 1 , xy 

Proof : 

Lemma (3.2) follows immediately from theorem (3.5) and definition 
(3.3). 



Lemma 3.3 : 

There can be no secondary path product terms generated from any 

path between x and £ in which branch b appears as an unbarred variable 

except path set P, 

1 »xy 

Proof : 

This lemma follows directly from theorem (3.5) and definition (3.4). 
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Theorem 3.6: 



Let path set, , be an A-edged path set between terminals x 

(e-A-1) 

and 2* This path will generate 2 secondary terms that are dupli- 



cates of those generated by P 



l,xy’ 



Proof : 

The number of possible unique combinations of the (e-A) number of 
barred variables has decreased by a factor of 2 because the 1-edged 

(e-1) 

path, P. , has previously formed 2 combinations of the same vari- 

l ,xy 

ables. Only those in which the variable b appears as a barred variable 

(e-A-1) (e-A) 

remain ungenerated. That number is 2 . There were 2 possible, 

(e-A-1) (e-A-1) 

but now only 2 , meaning that 2 were duplicates. 

Investigations have shown that the number of duplicates generated 
among all the path product terms are arranged in correspondence with the 
binary coefficients. That is, if there are eight duplicated generated, 
there will be a 1-3-3-1 distribution, based on the number of unbarred 
variables in the path. If there are D number of duplicates, then the 
distribution will be as given by 




D . 1 

j.' (D-j) 



j 0, 1,2, ... ,D 



The number of unique path product terms that are generated by 

each path set, P., in the 2-terminal reliability function, q , can 
i x,y 

be found by the following method. 

a) Write the path matrix, order (p x e) , where £ is the 
number of path sets between terminals x and y, and 

e^ number of edges in the graph. 

b) Rearrange the rows such that the shortest path is row #1, 
the next shortest is row #2, and so forth, until the last 
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row contains the longest path. 

c) Count the number of columns in the row in which there is 
no entry in any of the preceeding rows. For example, 
this number will be (e-k) for the first row. This num- 
ber may be zero. Call this number n^. 

d) Calculate 2°i. This result is the number of unique path 
product terms generated by the path set in the 2-terminal 
reliability function. 

The total number of terms in the Table of Combinations is the sum 
of the results of step (d) above. It is seen that if n^ is zero, that 
one term is still added to the Table of Combinations by the path, which 
agrees with Theorem (3.3). If there is only one path, then there will 
be 2 V 7 number of terms which agrees with Theorem (3.2). Further 
statements which tend to prove the validity of this method follow the 
presentation of this example. 



Example 3.4 : 




Path Matrix: 

edges : 
paths: i 

2 

3 




The paths between vertices A and B 



2 

0 

0 

1 



3 4 5 
0 0 1 
10 0 
0 10 
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After rearranging, the path matrix becomes: 



edges : 




1 


2 


3 


4 


5 


paths : 


2 


0 


0 


1 


0 


0 




1 


1 


0 


0 


0 


1 




3 


1 


1 


0 


1 


0 



Counting the number of columns in each row in which there is no entry 
in that column in any of the preceeding rows: 



v 



paths : 



2 

1 

3 



Calculating: 



paths : 



2 

1 

3 



4 

2 

0 

2 n i : 
16 
4 
1 



From definitions (3.1 - 5), there are established conditions which 

a collection of variables must meet in order to be an entry in the 

Table of Combinations. Looking at example (3.4), one can see that path 

4 

(2) will generate 2 unique terms. After they have been generated, path 

3 

(1) generates a collection of 2 number of possible terms. All these in 

which edge (3) appears as an unbarred variable will be duplicates of 

those generated by path (1). This reduces the number of terms by a 

2 

factor of 2, yielding the result: 2 . This can be understood from the 

concept of combinations; there were only 2 elements to be combined in 

2 

2 different ways vice what had appeared to be 3 elements; i.e., 2 vice 

2 3 . 
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The same argument holds for path (3). While there would appear 
to be 2 elements that can be varied, there are none which have not al- 
ready been varied. The fact that the primary path product term is unique 
itself is given in theorem (3.3). It is from this logic that the above 



stated procedure is derived. 

The knowledge of these numbers leads to a method to determine the 

value of the 2-terminal reliability function. With equation (3.5) in 

mind, a new equation can be written for each path set using the numbers, 

n found by the above described procedure. Let P be a k-edged path 
i x,y 

between terminals x and y in a graph with e number of edges. For this 
path, let n^ = N. Then one can write 








_(N-1) 

rr 



+ . . . 



(3.6) 




(N-l). 

r r 




where 



Q 



is the value of all the path product terms added to the 



Table of Combinations by the path, P^. 

Where the first factor is the product of the k number of branch 

reliabilities associated with the branches in the path, P ; the second 

x,y’ 

factor is the product of the (e-k-N) number of complements of the branch 
reliabilities of the branches which have previously appeared as unbarred 
variables in other path product terms; the third factor is the summa- 
tion of all the path product terms in which N number of variables, associ- 
ated with the branches not in the path, P , and not previously speci- 
fy 

fied as being unbarred in other path product terms, appear as barred and 
unbarred elements. 
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However, the quantity inside the brackets is equal to one. This 
statement is true because the quantity inside the brackets is an applica- 
tion of the binomial theorem for two conditions: the branch reliabilities 

being all equal; and the branch reliabilities being not all equal. 

Case 1) Let all the branch reliabilities be equal to r. 



Then the quantity inside the brackets of equation 



(3.6) becomes 



(r + r) 



N 



Since (r + r) = 1, then (r + r) 



N 



1 . 



Case 2) Let the reliability of branch, b^, be given by r^. 

Then the quantity inside the brackets of equation 



(3.6) may be rewritten: 



(r x + r 1 )(r 2 + r 2 ) . . . (r N + r N ) = (1) (!)...(!) = 1 



Therefore, equation (3.6) can be rewritten 




(3.7) 



In the event all branch reliabilities are equal 



k_(e-k-N) 
r r 



Q 



P, 



(3.8) 



k 
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4. Topological Analysis of the Net. 

4.1 Collection of All Possible Path Sets. 

Fundamental to the analysis of any communication net is the 
knowledge of all possible paths between the terminals x and £. Since no 
path may be ignored, even in relatively small, simple nets, the task of 
finding all paths is a laborous task. Linear graph theory provides an 
efficient method to accomplish this task. 

Consider a graph with e number of edges and n number of vertices, 
and define a circuit to be a closed edge train. Then, the theory of 
topology has shown there to be (e - n + 1) = w number of independnet cir 
cuits in the graph. These are chlled the fundamental circuit sets, and 
are denoted by B^. Any additional circuits that appear to exist in the 
graph can be shown to be algebraically derived from the set of indepen- 
dent circuits. 

Next there must be built a collection of sets, called set B, which 
includes the null set, the W number of sets of B^, and the collection of 
all possible combinations of circuit sets, which are generated by the 
ringsum operation of two or more circuit sets in B^. The ringsum opera- 
tion is defined as modulo 2 addition of the respective column entries of 

w 

the two sets. (See example, (4. 2) . The set B will contain 2 number of 
sets. 



A path set, P , is chosen, and expressed as a row matrix of e 
x,y 

columns. Set B is represented in matrix form of order (w x e) . The 

ringsum operation of P and each set in B produces a new set which is 

x,y 

written in a collection of sets called set BRP. This is also in matrix 
form of order (w x e). Linear graph theory shows that the ringsum of a 
path set and a circuit set will be either a new path set or the disjoint 



26 



union of a path set and a circuit set. Now each entry in BRP must be 
examined to retain only the pure path sets, which are written in set PFC, 
also in matrix form of order (p x e) , where £ is the number of paths 
between x and £. 

One process of examining set BRP is to look at each entry first to 
see whether the original path set, P , remains. If it does, and if 

x »y 

there are other branches present, then it is not a pure path set and must 
be discarded. If it does not, then it must be ascertained whether any 
of the circuits of set B are included. If so, then it is not a pure path 
set. Those that remain are the pure path sets, and more importantly, 
all the path sets between terminals x and A block diagram of this 
procedure is given in figure (4.2). 

The ringsum operation is built into the computer program by means 
of three FORTRAN IV instructions. The block diagram for this operation 
is shown in figure (4.3). Example (4.1) shows the ringsum operation for 
the four possible cases. The parameters must be logical elements; that 
is, either True, or False, with a 1 bit representing True, and a 0 bit 



False. 

Example 4.1: 


Case 1 


Case 2 


Case 3 


Case 


4 


A = 


0 


1 


i 


0 




B = 


1 


0 


i 


0 




C = 


1 


1 


0 


0 




The step-by-step 


computer program solution proceeds as 


follows : 


after step 1, C = 


0 


0 


0 


0 




after step 2, C = 


1 


1 


1 


0 




after step 3, C = 


1 


1 


0 


0 




The final results are identical. 
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Example 4.2: 




Figure 4.1 



Set B f : 

edge number 1 

Circuit number 1 1 

2 _0_ 

Circuit (3) = 

circuit (1) and (2) as follows 



Choose P 

i 

Set B: 




edge number 
set number 



1 

2 

3 

4 



1 

0 

1 

0 

1 



P AB ® B = BRP 



2 

1 

0 



1 

0 

1 



2 

0 

1 

0 

1 



e = 5 
m = 4 



w = 5- 4 + l = 2 



2 W = 2 2 - 4 



Circuit 

Circuit 



U) e 2 e 5 j 

(2) [<=3 e 4 ‘ 5 j 



3 4 

0 0 

1 1 

is derived 



5 

Order (w x e = 2 x 5) 
algebrically by ringsum of 




0 0 1 

111 

110 



3 4 

0 0 

0 0 

1 1 

1 1 



5 

0 

1 

1 

0 



Order (2 W x e = 4 x 5) 
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Set BRP 



edge number 
set number 



1 

2 

3 

4 



2 

0 

1 

0 

1 



3 

0 

0 

1 

1 



4 

0 

0 

1 

1 



Order (2 W x e = 4 x 5) 



Set PFC = 



edge number 
path number 



Paths are 



1 

2 

3 



1 

0 

0 






2 

0 

1 

1 



3 

0 

0 

1 



4 

0 

0 

1 



Order (p x e = 3 x 5) 



\ * ^®2 e 3 e 4^ 
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Block Diagram of Process to Find All Paths 

Figure 4.2 
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Block Diagram of Ringsum Operation in Computer Program 

Figure 4.3 
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4.2 Path Finding Program. 

A digital computer program can be devised to find all the paths 
between terminals x and The program discussed in Appendix (6) ac- 
complishes this task. (For further details, see that section). 

The generation of all possible combinations of the fundamental cir- 
cuit sets is an important step, and no one may be omitted. The method 
used by the program to find them all is the following: to count from 4 

to 7 in binary, one adds "100" to the binary numbers from 0 through 3. 
Similarly to count from 8 to 15, one adds "1000" to the binary numbers 
from 0 to 7. See example (4.2). This procedure guarantees that all 
possible combinations are obtained, because all numbers are obtained in 
the analogy. 

Example 4.3 : 

Decimal Binary 





0 = 


0 


0 


0 


0 




1 = 


0 


0 


0 


1 




2 - 


0 


0 


1 


0 




3 = 


0 


0 


1 


1 


to count from 4 to 7 


in binary, add 100 to 


those numbers 


preceeding 




4 = 


0 


1 


0 


0 




5 = 


0 


1 


0 


1 




6 = 


0 


1 


1 


0 




7 = 


0 


1 


1 


1 


Similarly, from 8 to 


15, add 1000 












8 = 


1 


0 


0 


0 




9 = 


1 


0 


0 


1 




10 = 


1 


0 


1 


0 




11 = 


1 


0 


1 


1 
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Decimal 



Binary 



12 = 



110 0 



13 = 



110 1 



14 = 



1110 



15 = 



1111 



4.3 Branch Loading. 

When a communication net is analyzed, each link is examined in- 
dividually. In addition to the reliability of the link to perform its 
function, there is needed a measure of the relative importance of the 
link in the various modes of operation of the complete net. For example, 

the importance of a link may be different, if complete communication is 

required, than if only K-terminal communication is required. One such 
measure of relative importance of a link is called the load of a branch. 

Fu has defined the load of a branch to be the total number trees 

that contain this branch, [3]. As shown before, complete communication 
is possible only when the net contains at least a tree in the operational 
state. Since each branch is contained in one or more trees of the net, 
the number of trees that contain this branch can be found, and this number 
used to indicate the relative importance of the branch in complete com- 
munication. This number can be normalized, if desired, since the total 
load of all branches can be quickly calculated by a formula shown by Fu, 
[3], as n 




1 i = (v - l)t 



i=l 



where 1. (i = 1, 2 
i 



m) is the load of branch 



in a network with v number of vertices, n branches, and t number 



of trees. 
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This method of branch loading determination can be extended to 
the k-terminal communication case by means of subtrees. However, branch 
loading can also be defined in another way. 

Definition 4.1 : Branch Loading . 

The load of a branch, b., is defined to be the total number of paths 

r 

(path sets) that contain this branch. 

A path set was defined to be a set of branches of a communication 
path between a pair of terminals. When complete communication is con- 
sidered, definition (4.1) makes use of all the paths between all pairs 
of terminals. This means that the numbers obtained for the various 
branch loadings will be larger than those obtained using the definition 
of Fu; however, the relative magnitudes will be the same. That is, the 
branch with the highest loading will be identified to be the same by 
either definition. The advantage of definition (4.1) lies in the simple 
way to evaluate rather than having to find the trees of the net. 

4.4 Path Length 

When a communication net is analyzed, all the paths between ter- 
minal pairs are found. The loading of the branches is measured. Now 
the paths themselves must be examined. The length of the communication 
path is of great significance. We know from experience that the more 
often data is handled, the greater the chance for error. This is con- 
firmed by this simple example. 

Example 4.4 : 

Let p be the reliability of all branches, b^, and p< 1. Then 

the reliability through one branch is p. The reliability of communica- 

2 

tion through two branches is ( p x p) = p ; through j number of branches: 
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p\ Since p < 1, then p j < p for all j > 1. 

This effect of path length appears directly in the 2-terminal re- 
liability function by means of the secondary path product terms. Section 
(4.6) discusses the effect of the secondary path product terms. See 
also section (5), Application of Ideas without the Aid of a Computer. 

4.5 Branch Reliability. 

Much work is currently being done by many investigators in the fields 
of Operations Research, Reliability, Non-Destructive Testing, and many 
others, to provide a high quality measurement of the reliability of a 
communication link. This work does not concern itself with the finding 
of this value of reliability. It must be pointed out, however, the signi- 
ficance of accuracy in this value and its effect on the interpretation 
of the results. 

Results have shown that considering a net with e = 10 branches, and 
assuming all branch reliabilities to be equal to .90, the fourth signi- 
ficant figure has risen to the importance level of being the deciding fac- 
tor in net configuration. 

4.6 Secondary Path Product Terms. 

The secondary path product terms of definition (3.4) are not of 

secondary importance in the determination of the reliability function. 

As can be seen from example (3.1), the majority of nonzero canonical 

terms are secondary path product terms. The exact number of these terms 

was given in theorem (3.1), namely, for an A-edged path, there will be 
( 6 -A) 

2 -1 number of secondary path product terms. 

The theorems in this section are subjected to the following condi- 
tions: A graph contains e number of edges, and n number of vertices. 

The branch reliability, p of branch b^, is greater than its ccxnplement. 
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p^, and is the same for all branches. 

Definition 4.2 : Value of a path product term . 

The value of a path product term is defined to be the product of 
the branch reliabilities of the unbarred variables and the complement 
of the branch reliabilities of the barred variables of the path pro- 
duct term. It is applicable to both primary and secondary path product 
terms. It always takes on a value less than one, since it is a product 
of numbers each less than one. 

Theorem 4.1 ; 

Any one of the secondary path product terms has a larger value than 
the primary path product term associated with it, for any A-edged path, 

Vy 

Proof: 

The value of a primary path product term is determined using de- 
finitions (3.3) and (4.2) to be 

A -(e-A) 

P • P 

By definition (3.4) a secondary path product term must have at least 

one more unbarred variable than the primary term. Therefore, we write 

(A + a) — (e-A-a) 
p • P 

where a is a positive integer, i.e., a ^ 1 

If p A . p (e - A) , ^ p (A + a > . 

Then 1 K. p a • p^ 

Since p < p, 

rn U A —(e-A) ^ (A + a) —(e-A-a) 

Thus p . p < p v ' p v . 
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Corollary 4.1 : 



The sum of the values of the primary and all the secondary path 
product terms of a 1-edged path is greater than that of any multi-edged 
path between the same terminals x and 



Proof: 

This Corollary follows directly from the theorem (4.1) by letting 
A be equal to 1. 

Theorem 4.2 : 

The sura of the values of all the primary and all the secondary path 
product terms of an A-edged path is greater than that of an (A+a) -edged 
path, where a is any positive integer; i.e., a ^ 1. 



Proof : 



Equation (3.5) can be written for an (A+a) -edged path. 
e-A-a 



T K(i)p (A+a+1) p (e - A - a - 1) = 
i = 0 




(A+a+1) - 
P P 



(e-A-a-1) 





(A+a)_(e-A-a) , 
P p + 



(A+a+2) -(e-A-a-2) . 

p p + • • ♦ (4. 1) 



+ 




p(e-l) 



«1 

P 



+ 




e 

P 



Comparing equation (4.1) to equation (3.5), we see that theorem (4.2) 
can be proved if it can be shown that, using definition (4.2), the value 
of equation (3.5) is greater than the value of equation (4.1). This can 
be done by showing that each coefficient of equation (4.1) is less than 
that of equation (3.5) for the corresponding power of £. 
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Consider the general term, i ** j . If 



e-A-a \ (e-j)-j 

e-A-a-j 1 p F 



< 



a-A^j I P <e ' j) P J 



(4.2) 



then 



e-A-a 

e-A-a-j 



< 



e-A 

e-A-j 



(4.3) 



then 



(e-A-a)! 



(e-A-a-j) ! (e-A-a (e-A-a-j)) 



< 



(4.4) 



(e-A)! 

(e-A-j) ! (e-A(e-A-j)) ! 



then 



(e-A-a) (e-A-a-1) (e-A-a-2) . . . (e-A-a-j 1) 

j.' 




(e-A) (e-A-1) (e-A-2) . . . (e-A-i 1) 

j! 



(4.5) 



There will be j. number of terms in the numberator of both sides of 
equation (4.5), each term of which is greater than one. 

Since, 



e-A-a 


< 


e-A 


e-A-a-1 


< 


e-A-1 


e-A-a-2 


< 


e-A-2 



e-A-a- j+1 ^ e-A-j+1 

Then equation (4.5) is true. From it equations (4.4), (4.3), and 
(4.2) follows, and the theorem is proved. 



38 



Example A. 5 : 



The following tables are the numerical values of the 2-terminal 
reliability function discussed in examples (3.1 and 2) Let the branch 
reliabilities be all equal and have the reliability 0.9. Columns 3,4, 
and 5 are the values when one branch is removed from the graph. Column 
6 is the value of the 2-terminal reliability function if the specified 
branch has the reduced reliability as shown. 





with all 


without 


without 


without 




with 




branches 


branch 


branch 


branch 




P =°. 5 




present 


#1 


#2 


#3 




X 


^x»y 


.981 


.81 


.90 


.90 




.905 




.981 


.90 


.81 


.90 




.905 


q 


.981 


.90 


.90 


.81 




.905 


M x,z 














The value of q^ 


.972 












The ideas of 


section (3.3) 


can be i 


used to find 


the value 


of 


the 


2-terminal reliability function 


in the i 


cases of examples (3.1 


and 2). 


Equation (3.8) is 


applicable in 


all the 


cases where 


the reliabilities 


of the branches are all equal; equation 


(3.7) where 


they are not 


all 



equal. 

In all of these examples there exists one 1-edged path and one 
2-edged path. This means, using equation (3.8), 
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y P = r = p .9 = .9 

| Q | P 2 = r 2 r = p 2 p = (.9) (.9) (. 1) = .081 

The value of q is .981 

x,y 

Observe that this is the same result obtained in example (4.5). 

Using equation (3.7) the case of example (3.1) where the reliability 
of branch one has the reduced reliability, p^ = 0.5. 



Q 

|Q 



r 2 r 3 r l 



P 2 P 3 f i = 



The value of q is 
x,y 



.5 

.405 

.905 



Observe that this is the same value found in example (4.5). It is signi- 
ficant to observe that in this example, if branch two had the reduced re- 
liability, the effect on the value of q would not be as much. Consider 

x >y 

example (3.1), let branch two assume the reduced branch reliability, p^ = 

.5, and calculate the value of the 2-terminal reliability function, q , 

* ,y 

using equation (3.7). 



Q 




r r r 

2 3 1 



The value of q 

x,y 

Observe that this value is 
branch one had the reduced 



■ Pj_ ■ .9 = .9 

= P 2 P 3 P 1 = <-5)(. 9 )(.l) = »045 

.945 

higher than that one found in the case where 
reliability: .945 vice .905. 
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5. Application of Ideas without the Aid of a Computer 

Some of the ideas presented in the preceeding sections can be ap- 
plied, without the benefit of a computer, in the analysis of a communica- 
tion net to give rough, but fast, measurements of the properties of a 
net. The following subsections discuss important properties of the 
communication net which must be examined in the event a change in the 
construction of the net is required. For example, the commander of 
several communication systems may order that one link be removed from one 
system and be given to a second system. The effect on the system caused 
by the removal of a link can be studied qualitatively without the aid of 
a computer. 

5.1 Branch loading. 

If the graph of a communication net is simple or if sufficient time 
is available, all the paths of the net can be identified, and the load- 
ing of each branch can be found. Immediately this indicates that the 
branch which is used most often should not be given up, because all the 
paths using that branch are also given up. On the other hand, the light- 
est loaded branch may not be sacrificed in every case without additional 
examination of the topology of the graph. That branch could be the only 
connection to a terminal, and the destruction of this branch would iso- 
late the terminal. However, in certain conceivable situations, the price 
of isolation of a terminal may be acceptable, for example, security 
reasons. The full effect of branch loading on the decision to remove a 
branch must be tempered by additional factors. 

5 . 2 Path Length . 

Theorem (4.2) proves that long paths are not as reliable as shorter 
ones. Therefore, one must place a figure of merit on the communication 
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path in relation to its length. The removal of either one of two branches 
may still enable communication between stations; however, in one case 
the remaining paths are longer than in the other. Then, the choice would 
be to remove the one that would leave the shorter path. This choice may 
be a poor one in light of the other factors which must be considered. 

5.3 Branch Reliability. 

Naturally the designer of a communication net will want to avoid the 
use of a branch with low reliability. And, if one link has to be removed, 
then that one of lowest reliability should be selected. However, this is 
not always true, especially when that branch is the only link to a cer- 
tain terminal. Again the price of isolation of that terminal may be ac- 
ceptable in a situation where security considerations are paramount. 
Likewise, the link of highest reliability may be chosen to be removed 
from the given net for use elsewhere, especially in cases where the 
branch loading is low and there exist other paths of communication of ac- 
c eptable reliability. 

5.4 Secondary Path Terms. 

The importance of the secondary path product terms has been stressed 
before. It is pointed out in this section that since they constitute 
the majority of terms in the 2-terminal reliability function, and they 
are inversely proportional in number to the length of path generating 
them, the loss of a short communication path is very costly in terms of 
reliability. 

The number of k-edged paths which will be lost when a given branch, 
b^, is removed can be predicted in certain cases. 
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Theorem 5.1: 



When branch, b^, is removed from a graph, one 1-edged path will be 
removed . 

Proof : 

There can be only one 1-edged path between terminals x and £ con- 
sisting of the branch b^ alone. Removal of the branch, b^, results in 
the removal of that path. 

Theorem 5.2 : 

Let and be the degree of terminals x and £ respectively. 

Let b^ be a branch connecting terminals x and Let c be the number of 
additional branches connecting terminals x and Then the number of 2- 
edged paths between any pair of terminals, but using branch b^, that will 
be removed if branch b^ is removed, is given by 

N 0 = (D + D - 2c -2) (5.1) 

2 x y 

Proof : 

The degree of a terminal is the number of edges incident to it. 

A 2-edged path using the branch, b^, must either start at terminal x 

and proceed through or start at terminal £ and proceed through x. 

In the former case, there are (D^-l) number of possible 2-edged paths 

minus the c number which would form a circuit, since they are connected 

to terminal x also. Similarly in the latter case, there are [(D^-lJ-c] 

number possible. The sum is N_: 

N 0 = D -1-c + D — 1 — c 
2 y x 

N = D + D - 2c -2 
2 y x 
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Example 5.1 




D = 5 
x 

D - 5 

y 

c = 1, namely 



Figure 5.1 

N 2 = 5 + 5 - 2(1) -2=6 



Namely, 

path sets 



^ b l b 3^ . [bjbj] 

and t 



• l b M • \wl 

l 




Rule 5.1 . 

A rule which establishes a series of operations that indicates the 
number of 3 -edged paths between any pair of terminals containing the 
branch, b^, that will be removed if the branch, b^, is removed. It is 
stated here without proof; however an example is given. 

Step #1 Starting at node x, proceed to the terminal of a branch 
incident to node x* 

Step #2 Call this node #1 having degree D^. Let c be the number 

of branches having terminals x and but not being branches, 

b^,. Let d^. be the number of branches having node ^ and 
node i as terminals. Then 

D £ + D y - 2 - c -2d i = N 3 (5.2) 

Step #3 Proceed to all other nodes of the branches incident to 

node x, except the node of course. Repeat the calcu- 

lation of equation (5.2), adding the result to the pre- 
vious value of N 3 . 
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Step #4 Proceed to the terminals of the branches incident 
to node ^ and whose terminal nodes are not those 
previously numbered or node x* Count the number of 
branches which are not incident to x or j or the nodes 
previously numbered in steps (2) or (3). Add this 
number of the previous value of N^. is the 

number of 3-edged paths that will be removed if 
branch, b., is removed. 

l 

Example 5.2: 



N 



3 




x 



Dj_ = 3 , d : - 0 

34-6-2-1-0 = 6 
D 2 = 3 , d 2 - 1 

3+6-2-1- 2(1) = 4 



16 

11 



N 



3 



N 3 = 6 + 4 = 10 
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1 




d 



3 



2 + 6 - 2 - 1 - 2(1) = 3 



N. 



3 



N, 



3 



10 + 3 - 13 



D. = 3 
4 

3+6 



• V 1 

2 - 1 - 2(1) = 4 



N. 



3 



N 3 = 13 + 4 * 



Next node 5 where 2 branches qualify 



N 3 = 17 + 2 = 19 



ANS : 19 



5.5 Number of Paths 

Theorem (3.3) proves that the primary path term will not be a 
duplicate of any entry in the Table of Combinations regardless of the 
number of paths that exist between terminals x and £. Therefore, each 
path will contribute some value, in the sense of definition (4.2), to 
the 2- terminal reliability function, , and its removal will deminish 

this value. 

An examination of the paths between terminals x and £ alone is not 
sufficient when n-terminal reliability is concerned. What may appear to 
be an expendable path between terminals x and £, may be a vital path 
between two other terminals. 

While a large number of paths between the terminals of considera- 
tion indicates that the loss of one path would not cause isolation, the 
length of the paths may be such as to result in low reliability. All 
the factors mentioned previously continue to apply even when there are a 
large number of paths. This is also true when there are a small number 
of paths. 
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5.6 Addition of a branch. 

The addition of one branch to a graph is equivalent to the addi- 
tion of one communication link to a net. This can be viewed as a de- 
sign problem in which there are eight cases to be studied: 

What should the second branch terminal be, if the given objective 
is to: 



maximize q 



x,y 



maximize q 



maximize 





the first branch terminal 
required to be: 


is 


(1) 


terminal 


X 




(2) 


terminal 


a, not x or 2 




(3) 


optional 






(4) 


terminal 


a in set k 




(5) 


terminal 


a not in set k 




(6) 


optional 






(7) 


terminal 


X 




(8) 


optional 






many more 


if the reliability of the 


branch 



to be added and the reliabilities of the other branches are permitted to 
change. The following discussion of the above cases assumes that all 
branch reliabilities are the same, and equal p. 



Case (1 ) : 

Corollary (4.1) specifies terminal 2 as the second branch terminal. 
Case (2 ) : 

Find all the paths between terminals x and a, and 2 and a. Choose 
terminal x or 2 depending on which has the fewer paths to terminal a. 
This will provide the most paths, with the fewest branches, between x 
and 2* 
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Case (3 ) : 

Corollary (4.1) specifies the selection of terminals x and 
Case (4 ) : 

Find all the paths between all the pairs of terminals of set k. If 
any one pair has an unusually small number, examine the topology of the 
graph in light of sections (5.1-4). Base the selection on the unusual 
topology. Do the same thing in the event one pair has an unusually large 
number of paths between them. Otherwise, locate the terminal that has 
the most number of paths to terminal a. Select this as the second ter- 
minal for the branch to be added. This choice will add one 1-edged path, 
but also the maximum number of 2- and 3-edged paths between the other 
terminals in the set k. This will increase the value of the most. 

Case (5 ) : 

Find all the paths between all the pairs of terminals of set k and 
terminal a. Examine the topology of the graph in light of sections (5.1- 
5). Especially, look for vertices of degree one or two. Unusual con- 
figurations must be solved in this manner. Otherwise, select as the 
second terminal that one in set k that has the largest number of paths to 
terminal a. See Case (4) . 

Case (6) : 

Find all the paths between all the terminal pairs in set k. Choose 
as the terminals of the branch to be added those between which there al- 
ready exist the largest number of paths. This will not be true if an 
examination of the topology reveals an unusual nature to the net. See 
Case (4). 

Case (7) : 

Find all the paths between all the terminals of the graph and ter- 
minal x. Examine the topology of the graph. If nothing extradorinary 
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appears, choose as the second terminal for the branch to be added that 
terminal which already has the most paths to terminal x. 

Case (8 ) ; 

Find all the paths between all the terminal pairs of the graph. 

Study the topology of the net. Choose the terminal pair between which 
there already exist the largest number of paths. 

If there exist a large number of paths between two terminals of a 
graph, the paths are usually of long length. This means that they rely 
on many other branches, and so are of low reliability. The adding of a 
1-edged path between these terminals immediately increases the 2-termin- 
al reliability between them, and also provides an increase in the number 
of 2-edged paths between other pairs of terminals. In fact, it increases 
the number of all k-edged paths. Since there were all those paths in 
existence, the accessibility of the new branch by the other terminals 
produces a maximum number of paths between all terminals. 

Lemma (3.1) shows that an isolated terminal, I, will cause both 
the 2- and the n-terminal reliability function to be zero. When an 
isolated terminal exists, clearly the topology of the graph is unusual. 
The isolated terminal must be selected as the second terminal for the 
branch to be added in each case that the terminal I is in the set k. 

This statement covers all the cases since the set k includes all situa- 
tions from 2 to n number of terminals. 

Example 5.3 : 

Three illustrations of a 2-terminal reliability function are 

discussed for the graph of figure (5.3)? q , q , q . All the 

ar b ,g c , g 

branch reliabilities are equal, p = .9. 
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i) 






with all 




without 


without 


without 




branches 




branches branches 


branches 




present 




#3 


#6 


#9 


q a,f 


.9859 




.9685 


.9755 


.9750 


q b,g 


.9940 




.9621 


.9922 


.9897 


q 


.9963 




.9854 


.9937 


.9830 


c,g 












Branch Loading: 












Branch 


number: 1 


2 


3 4 5 


6 7 8 < 


) 10 


q f Branch 

a , r 


load : 6 


4 


4 4 4 


3 5 3 4 4 10 paths 


q a,g 


1 


3 


4 3 3 


3 3 4 : 


) 1 8 paths 


q. _ 


3 


4 


3 4 3 


433 : 


} 3 9 paths 


Total loading: 10 


11 


11 11 10 


10 11 10 10 8 



In the illustrations given, branch #3 is not the one to be removed. 
However, branch #10 should not be removed either because this would cause 
node A to rely completely on branch #1 for communication. The low level 
of branch loading for branch #10 indicated that it is ineffectively 
connected, not that it can be sacrificed. 
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6. Digital Computer Program. 

6.1 Program Specifications. 

A digital computer program has been written in the FORTRAN IV 
Language suitable for use on the IBM System 360 computer complex. The 
primary purpose of the program is to compute the 2-terminal reliability 
function, a . This program is included as Appendix I . The func- 

^,y 

tional block diagram, figures (6.1) > (6.2) and (6.3) is fully explained 
in section (6.2). 

The capabilities of the program go beyond the computation of the 2- 

terminal reliability function, q . After finding q , the program 

^»y x,y 

can calculate its numerical value; then, it can recalculate it under 
the condition that a branch of the graph has been removed, and that a 
branch has a different value of reliability. (See options. Section 
6.1.1.) After the second and any succeeding runs, (a run is defined as 
the computation of one 2-terminal reliability function, )q^, where k is 
the set of terminal pairs of the preceeding runs, can be computed, and 
its numerical value found. In finding q , all path sets between 

x »y 

terminal pairs are found; this information is given as a preliminary out- 
put. 

6.1.1 Input 

The input required by the computer program consists of three types: 

(1) The number of edges, and the number of vertices of the graph 
are the first data required. The number of fundamental circuits is cal- 
culated from this data using the relation: 

(e - n + 1) = w 

Also on the first data card, the option SEL1 is provided which, if used, 
enables the branch reliabilities to be input. This will allow the 
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calculation of the numerical value of q , and later, the numerical 

*,y 

value of q^» if requested. 

(2) The w number of fundamental circuits are the next input 
data. If the branch reliabilities are to be used, they are read in next. 

(3) The data required by each run is read next. This data 
consists of three main parts: 

a) the terminal vertices (numbered) ; 

b) one path between these vertices; 

c) option selections: 

1) SEL - additional runs will follow 

2) SEL2 - find q fc 

3) SEL4 - calculate the numerical value 

4) SEL5 - of q with the specified branch 

*,y 

5) SEL6 - having been removed from the graph 

6) SEL7 - calculate q^ ^ using the specified value 

of branch reliability for the specified 
branch 

Precise details of the format for the data cards are provided in sec- 
tion (6.3). 

6.1.2 Output. 

The output provided by the computer program is fully labelled for 
easy reading. All the input data is printed with the generated data to 
provide a complete record of the run. The output of generated data in- 
clude : 

1) Table of Combinations; 

2) path matrix; 

3) branch loading information; 
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4) number of k-edged paths; 

5) numerical values of q and q, . 

At the conclusion of the computer program, a summary of the branch load- 
ing and the number of k-edged paths is printed. A message explaining 
the error stop is printed as necessary. 

6.1.3 Theoretical Limitations. 

In its present form, the computer program is limited to a graph 
containing no more than 15 edges, or 10 nodes, or 10 fundamental cir- 
cuits. There must be at least one circuit in the graph. Other limita- 
tions include a maximum of 4000 entries in the Table of Combinations for 
all runs, and a maximum of 210 runs. The reason for these limits is the 
storage requirement, especially that for the Table of Combinations. Dif- 
ferent programming techniques might be able to increase these limits 
slightly. The fact remains that the powers of 2 increase very rapidly, 
regardless of the programming technique. 

6.2 Block Diagram. 

The functional block diagram, figures (6.1,2 and 3), is presented 
to give a better description of the operation of the computer program. 

The blocks are identified in the listing of the computer program. Ap- 
pendix I. 
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Figure 6.1 
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12 




Figure 6.2 
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Figure 6.3 
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Description of Block Diagram 



Block 0: 


Control cards and dimension and format statements are required 



to prepare the computer for running. 



Block 1: 


Input data are read including number of branches and nodes, 



and the fundamental circuits. 



Block 2: 


Graph size is tested to be within program limits. 


Block 3: 


Input data is converted to logical elements to facilitate 



faster calculation. 



Block 4: 


If branch reliabilities are to be used, they are read in 



now into DREL matrix. 



Block 5: 


Preliminary printouts are made for more complete data refer- 



ence for the final output. 



Block 6: 


Values required for first run of the program are set. 


Block 7: 


Values required for succeeding runs are set. 


Block 8: 


Data required for the run is read including terminals and 



path . 



Block 9: 


Path data is converted to logical elements, and a storage 




area prepared for later use. 



Block 10 : All possible combinations of the given circuits are generated 
and stored in matrix BFC. 

Block 11: The ringsum operation is performed on the given path and 





each entry in matrix BFC, and the result is written in BRP. 


Block 12: 


The original path is written into the path matrix, PFC. 


Block 13: 


The primary and all the secondary path product terms are 



generated into OPC . 



Block 14: 


The path product terms are printed out. 


Block 15: 


Matrix BRP is searched for pure path sets. 
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Block 16: 



Block 17 : 



Block 18 : 

Block 19 : 
Block 20: 



Block 21 : 

Block 22 : 
Block 23 : 

Block 24 : 

Block 25 : 
Block 26: 



Block 27: 



The primary and secondary path product terms are generated 
for the newly formed path set, and written into QTEMP. 

Each entry in QTEMP is compared to all entries in QPC. 

If a duplicate is not found, the set from QTEMP is written 
in QPC . 

Those sets added to QPC are printed out. The number of 
entries in QPC and the number in PFC are printed. 

The path matrix, PFC , is printed out. 

The various OPTIONS are examined. If Q. . or Q are to be 

i » J n 

found, transfer is made; if another run it to be made, trans- 
fer is accomplished. 

Before stopping, a summary is printed out including branch 
loading, number of paths between each node, and path length. 
Error stops are provided. 

Generate matrix Q from path product terms common to the first 



2 q, matrices. 
i,v 

Generate matrix Q from the previous Q and the newly formed 
m n 

q . for this run. 

H i,j 

Print out matrix Q . 

n 

If the branch reliabilities are provided, determine the value 



of Q^. Transfer control to block 20 to check options. 

If the branch reliabilities are provided, determine the value 
of q^ j. Examine the options SEL4, SEL5, SEL6, and determine 
the value of ^ without the specified branch J. Examine 
option SEL7, and compute the value of q^ ^ with the revised 
value of branch reliability. Transfer control to block 20 to 
check other options. 
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6.3 Data Card Configuration. 



The following is a description of the data cards required for the 
running of the computer program as described in section (6) . 



Data Card 

Type Number : 
1 



Columns : 



1 - 3 



4 - 6 



7 - 9 



10 - 72 



Description 

and specifications : 

Number of edges, right justified; 
max: 15 

Number of nodes, right justified; 
max: 10 

If 0 or Blank: no branch reliabilities 

are to be used; if nonzero: branch 
reliabilities are to be read and used. 
Disregarded. 



2 



3 



1 - 15 



15 - 72 
1 - 5 

6 - 13 

14 - 21 

22 - 29 

30 - 37 

38 - 45 

46 - 72 



Fundamental circuit set: 

Up to 15 digits, 0 or 1, left 

justified, representing the edges in 

the circuit; 

1 if edge present 

0 if edge not present 

Disregarded. 

Blank 

Decimal point and 7 digits to represent 
the branch reliability of the succes- 
sive branches 1 , 2, . . . , e 
One field for each edge. Up to five 
edges per card. 

Disregarded. 
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Type Number: Columns: 


Description 

and Specifications: 


4 1-3 


Identifying number of starting node 


4-6 


Identifying number of terminal node 


7-9 


if 0, no more cards will be read; 
if not 0, additional cards read. 


10 - 12 


if nonzero, will be computed using 

the previously found q. from all 

1 » J s 

previous runs as the set k; if 0, 
then this is not done. 


13 - 27 


left justified, up to 15 digits, 

0 to 1, representing the edges in a 
path between the subject terminals; 

0 if edge not present 

1 if edge present 

only e number of digits processed 


28 - 30 


an integer, less than or equal to e, 

representing an edge that has been 

removed from the graph. A value for 

q. .is found assuming that edge is 
1 9 J 

removed from the graph. 


31 - 33 


another edge number processed as that 
on in the preceeding field. 


34 - 37 


Same as above 


38 - 40 


An edge number whose value of branch 
reliability it is desired to change 
temporarily, and recompute q^ ^ with 
the new value of reliability. 
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Type Number : Columns : 

41 - 48 



Description 

and Specifications : 

A decimal point and 7 digits represent- 
the branch reliability of the edge 
specified in the preceeding field. 
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7 . 



Conclusions. 



This work shows that the reliability function as defined is an ef- 
fective instrument in the analysis and design of communication nets. The 
theorems that are developed can be applied in many cases without the use 
of a digital computer to provide direction in the modification of a net. 
The computer can provide fast data to quantify the effect of any change 
to a communication net with respect to k-terminal reliability, as the 
results of the computer study indicate. 

The computer program can be modified by specialization to facilitate 
sensitivity studies for branch removal or branch addition. In so doing, 
the limitations on graph construction can be modified. Moreover, equa- 
tion (3.7) can be used in a computer program to study the 2-terminal 
reliability function alone. 
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APPENDIX I 



The digital computer program, as described in section (6) is 
presented as Appendix I. The programming language is FORTRAN IV and 
the program is suitable for use on the IBM System 360 computer complex. 



0001 

0002 



0003 

0004 



0005 

0006 

0007 

0008 

0009 

0010 

0011 

0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 



C 

C BLOCK NUMBER 0 

C 

IMPLICIT L0GK.AL*1 ( A-C.O-O) , 

1 INTEGER*4(E-N,P-Z) 

DIMENSION XBF( 10,15) , BC.MI 10,15) , 

1BFCI 1024,15) ,BRP( 1024, 15) ,PFC( 1024, 15) , 
1XPFCI (210,15) , XY ( 2 10, 15) ,ENO( 210, 15), 
1GN0 (210) ,QNO (4C0,1 5) ,QN1 (400,15) , 

10TSI 15) , OPC ( 4000 ,15) ,QTEMP(400, 15) , 
llJNOl 210) , URUN (210) ,XPTH{ 15.) , 

1PTS! 15) ,ST(210),SP<210) ,DREL( 15) 



45 FORMAT!/, 26H THE VALUE OF Q-SUB-(I,J) 

115H WITHOUT BRANCH 113) 

46 FORMAT!/, 26H THE VALUE OF Q-SUB-(!,J) 

1 1 4H WITH BRANCH 113,/, 

125H HAVING RELIABILITY VALUE1F8.6) 

76 FORMAT! /,6H S T 0 P 1 I 5 ) 

77 FORMAT!/, 22H TOTAL BRANCH LOADING) 

78 FORMAT!/, 25H IN PAIRS, THE NUMBER OF 
1 33 H K-EDGEI) PATHS AMD BRANCH LOADING) 

79 FORMAT ( IH1 , 25H THE TOTAL NUMBER OF PAT 
120HHS FOUND IN EACH RUN) 

80 FORMAT! 313) 

81 FORMAT! 1511) 

82 FORMAT! 216, 151 1,112) 

83 FORMAT! 5X, 1513 ) 

84 FORMAT! 1H1,20H THF GIVEN PATH WAS) 

85 F0RMATI5X,/) 

86 FORMAT! 1H1.1CH THERE ARE 115, 

1 1 0 H LINKS AND 115, 6H NODES/) 

37 FORMAT!/, 25H f DGE 1 2 3 4 5 6 

125H7 8 9 IP 11 12 13 14 15/) 

88 FORMAT!/, 29H THE GIVEN CIRCUIT MATRIX IS) 
39 FORMAT ( / , 2CH THE PATH MATRIX IS) 

90 FORMAT ( / , 5H R=1I3,5H U = 1 13) 

91 FOR MAI ( / , 20H MATRIX OPC FOLLOWS) 

92 FORMAT (413,1 511,413, 1FB. 6) 

93 FORMAT! 5X, 5F8. 6) 

94 FORMAT ( / , 30H THF BRANCH RELIABILITIES ARE) 

95 FORMAT!/, 13H C-$UB-I,J =1F10.7) 

96 FORMAT!/, 11H C-SU6-K =1F10.7) 

97 FORMAT!/, 19H MATRIX OK FOLLOWS) 

98 FORMAT! IH1 ,8H RUN NO . 1 I 6 , 

16H FR0M1 15, 4H TO 1 I 5 , / ) 

99 FORMAT!/, 10H THERE ARE1I5,6H PATHS) 
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0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

005? 

0053 

0054 

0055 

0056 

0057 

0058 

0059 

0060 

0061 

0062 

0063 

0064 

0065 

0066 



non ooo non 



BLOCK NUMBER 1,2, 3, 4, AND 5 



RFAD<5,80) E » N » S F. L 1 
WRITE (6,86) E,N 
W=F-N+ 1 
I F ( E.GT.15) 

IF( E.LE, 0) 

IF(N.GT.IO) 

I F ( N • LE • 0) 

IF(W.GT.tO) 

IF(W.LE. 0) 

00 71 1=1, W 

71 READ( 5,81) ( X 8 F ( I »J) ,J=1 »E) 
00 75 1=1, W 

00 70 J=1,E 
IF ( XBF ( I , J ) . FQ . 0 ) 

1 F ( X BF ( I , J ) .FO. 1 ) 



GO TO 3 
GO TO 3 
GO TO 3 
GO TO 3 
GO TO 3 
GO TO 3 



BCM( I , J) = . FALSE, 
B C M ( I , J) = .TRUE. 



70 CONTINUE 
75 CONTINUE 
WRITE(6,88) 

00 72 1=1, W 

72 WRITE (6, 83) ( XGF ( I , J) , J= L, F ) 

WRITE ( 6,85) 

IF(SELl.EQ.O) GO TO 10 
RFAD(5,93) (DRFL(J) » J= 1 , E ) 

WR I T E ( 6 , 94 ) 

WRITE(6,93)(0RFL( J),J=1,F) 

BLOCK NUMBER 6 AND 7 

10 T I M S = 0 -* 

OPT = .FALSE. 

OTAG =. FALSE. 

QTAG1=. FALSE. 

QTAG2=. FALSE. 

GR = 0 
R = 0 
U = 0 

1 R=R+1 
URAN=U+ 1 
U=0 
UV = C 

BLOCK NUMBER 8 AND 9 

RFAD( 5,92)ST(R) ,SP(R) , SEL , SEL2 , XPTH , 
1 SEL4, SEL5»SEL6,SEL7,0CHG 



0067 

0068 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 

0078 

0079 

0080 

0081 

0082 

0083 

0084 

0085 

0086 

0087 

0C88 

0089 

0090 

0091 

0092 

0093 

0094 

0095 

0096 

0097 

0098 

0099 

0100 

0101 

0102 

0103 

0104 

0105 

0106 



on ooo ooo 
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WRITE(6 f 84) 

WRITE(6,87) 

WR ITET6 » 83 1’ XPTH 
WRITE (6,85) 

00 54 I = 1»E 
IF (XPTH(I).EQ.O) 
IF ( XPTH( I ) . E0. 1 ) 
XY ( - 

ENG ( R » J ) =0 

BLOCK NUMBER 10 
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IF (QPT)GQ TO 51 
OPT = .TRUE. 

DO 55 1=1, E 
8FCUt !•> = . FALSE. 
BFC(2, I )=BCM< 1, I ) 
K = 2 
G=3 



P T S ( i )=. FALSE. 
PTS< I )= . TRUE. 



IF(W.EQ.l) GO TO 51 

DO 58 I=2,W 

DO 57 J=1,K 

DO 56 F = 1,E 

A = BFC ( J » F )- 

8=BCM( I , F ) 

C =. FALSE. 

IF (A.OR.B) ~C- . TR UE . 
IF (A.AND.B) C=. FALSE. 
BFC(GtF) = C 

56 CONTINUE 
G=G+ 1 

57 CONTINUE — 

58 K=G- 1 

BLOCK NUMBER 11 



51 T=GR+ 1 

DO 38 I = 1 1 K 
DO 37 J=1,E 
A=BFCII,J) 

52 8=PTS ( J ) 

C =. FALSE. 

IF (A.OR.B) C= . TRUE . 

IF (A.AND.B) C=. FALSE. 

37 3RP ( I t J ) =C 

38 CONTINUE 

BLOCK NUMBER 12t 13t AND 14 
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0107 

0108 

0109 

0110 

0111 

0112 

0113 

01 14 

0115 

01 lb 

0117 

0118 

0119 

0120 

0121 

0122 

0123 

0124 

0125 

0126 

0127 

0128 

0129 

0130 

C 1 3 1 

0 132 

0133 

0134 

0135 

0136 

0137 

0138 

0139 

0140 

0141 

014 ? 

0143 

0144 

0145 

0146 

0147 

0148 

0149 

0150 



ooo 



U=UR AN 
UV=U 

90 103 J = 1,E 
QTS( J)=PTS( J) 

JPC(U,J) = PTS( J) 

193 PFC ( T , J ) = PTS ( J ) 

00 115 J = 1,E 

1 F ( QTS ( J ).) GO TO 115 
0T S ( J ) = . TRUE . 

USTOP=UV 

DO 102 V=U » US TCP . 

UV=UV+1 
00 101 JJ = 1 » E 
C=. FALSE. , 

A =O p C ( V 1 0 J > 

3 = QT S I J J ) 

IF(A.OR.B) C=.TRUE. 

QPC < UV » J J ) =C 
1C1 CONTINUE 
102 CONTINUE 

0TS( J )=. FALSE. 

115 CONTINUE 
IJ=UV 

WRITE<6,91) 

WR I TE ( 6 » 87 ) 

H=0 

00 225 HO=UP AN » UV 
H = H+ 1 

DO 224 J = 1 , E 
XPFCIIH, J)=l 
IFIOPCIHOtJ)) GO TO 224 
XPFCI ( H , J ) =0 
2 ?4 CONTINUE 

WRITE (6, 83) ( XPFC I < H , J ) , J= l , E ) 

225 CONTINUE 

BLOCK NUMBER 15 

DO 228 1=1, K 
AD=. FALSE. 

'90 226 J=1 ,E 
I F ( PT S ( J ) ) GO TO 230 

226 CONTINUE 

IF(AD) GO TO 228 
GO TO 236 

2 30 I F ( 8RP ( I , J ) ) GO TO 234 
GO TO 236 
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\ 



0151 

0152 

0153 

0154 

0155 

0156 

0157 

0158 

0159 

0160 
0161 
0162 

0163 

0164 

0165 



0166 

0167 

0168 

0169 

0170 

0171 

0172 

0173 

0174 

0175 

0176 

0177 

0178 

0179 

0180 
0181 
0182 

0183 

0184 

0185 

0186 

0187 

0188 

0189 

0190 

0191 

0192 

0193 

0194 



234 AD= . TRUE . 

GO TO 226 

C CIRCUIT EXISTENCE TEST 
236 AD=. FALSE. 

00 299 L= 1 » K 
AD=. FALSE. 

DO 295 J = 1,E 
IF(BFC(L,J> > GO TO 298 
295 CONTINUE 

IF (AD) GO TO 228 
GO TO 300 
238 AD= . T RUE . 

GO TO 295 

298 IF ( BRP ( I, J) ) GO TO 238 
300 CONTINUE 

299 CONTINUE - 
C 

C BLOCK NUMBER 16, 17 AND 18 

C 

T = T+ 1 
V= 1 

DO 104 J = 1 , E 
QTS( J)=BRP( I ,J) 

0 T E M P ( 1 , J ) =BRP ( I,J ) 

104 PFC(T,J) = BRP( I ,J) 

DO 108 J= 1 , E 
IF(QTSIJ) ) GO TO 108 
QTS(J)=.TRUE. 

STEP = V 

DO 106 IV=1,STFP 
V=V+1 

00 105 J J= 1 , E 
C=. FALSE. 

A = QTEMP ( IV, JJ) 

B = QTS ( J J ) 

IF(A.OR.B) C= . TRUE . 

QTEMP ( V , J J ) =C 

105 CONTINUE 

106 CONTINUE 

QTS( J>=. FALSE. 

108 CONTINUE 
KT R=0 

00 113 H= 1 , V 
DO 111 I V=UR AN , U 
DO 110 J = 1 , E 
A=QTEMP ( H , J ) 

B = QPC( IV, J ) 

IF(A.ANO.B) GO TO 110 
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019 5 

0196 

0197 

0198 

0199 

0200 

0201 

0202 

0203 

0204 

0205 

0206 

0207 

0208 

0 2C9 

0210 

0211 

0212 

0 213 

0214 

0215 

0216 

0217 

0218 

0219 

022C 

0221 

0222 

0223 

0224 

0225 

0226 

0227 

0228 

0229 

0230 

0 231 

0232 

0233 

0234 

0235 

0236 

0237 



non 



GO TO 111 



I F ( A , OR • 8 ) 

110 CONTINUE 
GO TO 113 

111 CONTINUE 
KTR=KTR+ 1 
!JV=UV + 1 

00 112 J= 1 » E 
XPFCI ( 1,J)=1 

IF(QTEMP(H,J) ) GO TO 112 
XPFCI < 1, J)=0 

112 9PC(UV,J>=QTEMP(H, J) 

WRITE<6,83) < XPFCI ( 1, J) ,J=1,E) 

113 CONTINUE 
U=UV 

1 F { U .GT . 4000 ) GO TO 5 
228 CONTINUE 

TNG=T-GR 
WPITE(6* 85) 

WRITE (6,99) TNG 
WR ITE (6,89) 

WP I T E ( 6 » 87 ) 

8 LOCK NUMBER 19 

TSTR = GR + 1 
H = 0 

DO 127 I =TSTP,T 
W = H+ 1 
X-0 

00 126 J = 1,E 
XPFCI (H,J) =1 

1 F ( P FC < I,J) ) GO TO 128 
X P FC I (H, J ) = 0 

126 CONTINUE 

WRITE (6,83) ( XPFC I (H, J ) , J=1 ,E) 

XV ( R , X ) =XY( R , X ) +1 

127 CONTINUE 
WRITF(6,78) 

WR I TE ( 6 , 87 ) 

WRITE (6, 83) (XY(R t X)tX=l,E) 
WRITE(6,83) (ENP(R, J) , J=1,E) 
WRITE (6,85) 

UR UN ( R ) =U 
UNO( R )=U-URAN+1 
GNO ( R ) =TNO 
GR = T 

C BLOCK NUMBER 20, 21 AND 22 
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\ 



c 



0238 






IF( SEL1.NE.0) GO TO 15 


0239 




119 


IF(OTAG) GO TO 120 


0240 






OT AG= .TRUE. 


0241 




121 


IF (SEL.EQ.O) GO TO 53 


0242 






GO TO 1 


0243 




128 


ENO(R,J)=ENO(R, J)+l 


0244 






X = X+ 1 — - 


0245 






GO TO 126 


0246 




53 


WR I TE ( 6 * 79 ) 


0247 






WR ITE ( 6 , 83 ) t GNC ( I) ,1 = 1, R) 


0248 






WRITE(6,78) 


0249 






DO 6 J=1 »E 


0250 




6 


UNO ( J ) =0 


0251 






00 2 I = 1 , R 


0252 






WR I TE < 6 , 98 ) I ,ST< I >, SP< I ) 


0253 






DO 7 J= 1 » E 


0254 




7 


UNO (J)=UNO(J)+ENO( I ,J) 


0255 






WRITE (6,83) ( XY( I , J) , J = 1 ,E ) 


0256 






WRITE (6, 83) ( ENO ( I, J) ,J=1,E) 


0257 




2 


WR I TE ( 6 , 85 ) 


0258 






WRITE(6,77) 


0259 






WR I TE ( 6 , 83 ) ( UNC ( J) ,J = l,E) 


0260 




9 


STOP 


0261 




3 


RE STP=99 


0262 






GO TO 8 


0263 




4 


RE ST P=999 


0264 






GO TO 8 


0265 




S 


RESTP=8S8 


0266 




8 


WRITE(6,76) RESTP 


0267 


r 




GO TO 9 




L 

c 

r 




BLOCK NUMBER 23 


0268 




120 


IF ( SEL2.E0.0) GO TO 121 


0269 






IF(QTAGl) GO TO 15C 


0270 






QTAG1=.TRUE. 


0271 






Z = i 


0272 






Z A 1=UN0 ( 1 ) • 


0273 






ZA2=ZA1+1 


0274 






ZA 3=Z A 1+lJNO ( 2 ) 


0275 






00 125 YZ=1»ZA1 


C 2 76 






00 124 YY=ZA2 , ZA3 


0277 






00 123 J = 1 » F 


0278 






A=QPC ( YZ , J ) 


0279 






B=QPC (YY, J ) 


02 80 






I F ( A . AND . B ) GO TO 122 


0281 






IF (A.OR.B) GO TO 124 
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0282 

0283 

0284 

0285 

0286 

0287 

0288 

0289 

0290 

0291 

0292 

0293 

0294 

0295 

0296 

0297 

0298 

0299 

0300 

0301 

0302 

0303 

0304 

0305 

0306 

0307 

0308 

0309 

0310 

0311 

0312 

0313 

0314 

0315 

0316 

0317 

0318 

0319 

0320 

0321 

0322 

0323 

0324 

0325 

0326 



122 QN0(Z,J)=A 

123 CONTINUE 
Z = Z+l 

124 CONTINUE 

125 CONTINUE 
YR = 2 

GO TO 161 
C 

C 3 LOCK NUMBER ?4 AND 25 

C 

150 ZSTR=URUN< YR)+1 
YR=YR+ 1 
ZMX=Z-1 
ZAR=URUN ( YR ) 

Z = 1 

00 160 YZ=ZSTR t Z AR 
DO 159 YY= 1 » Z MX 

00 158 J = 1 » E 
A=QPC<YZ, J) 

I F ( 0TAG2 I GO TO 154 
B = QNO ( YY * J ) 

151 I F ( A . AND . B ) GO TO 152 
IF(A.OR.B) GO TO 159 

152 IFIQTAG2) GO TO 156 
ON 1 ( Z » J ) =A 

GO TO 158 
154 3=QN1(YY,J) 

30 TO 151 
156 ONO ( Z » J ) = A 
3C TO 158 

158 CONTINUE 

1 = 1 * 1 

159 CONTINUE 

160 CONTINUE 
OTAG2=.NOT.QTAG2 

161 IF(YR.EQ.R) GO TO 162 
30 TO 150 

162 7^X=Z- 1 
WRITE (6, 97) 

DO 27 YD = 1 » ZMX 

DO 26 J= 1 . E 

XPFCI I 1, J) = l 

I F ( OT AG2 ) GO TO 24 

I F ( QNO ( YD t J ) ) GO TO 26 

25 XPFC I ( 1 1 J ) =0 

26 CONTINUE 

WRITE (6, 83) (XPFC I ( 1, J) , J*1,EI 

27 CONTINUE 
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0327 






GO TO 164 


0328 




24 


IF ( QN 1 ( YD » J ) ) GO TO 26 


0329 


c 




GO TO 25 




c 

r 




BLOCK NUMBER 26 


0330 




164 


IF< SEL2.NE.2) GO TO 121 


0331 






ZMX=Z-1 * 


0332 






DSUM=0. 0 


0333 






DO 168 YD= 1 1 ZMX 


0334 






DTEMP=1'.0 *' 


0335 






00 167 J=1,E 


0336 






I F ( 0TAG2 J GO TO 165 


0337 






IF ( QNO ( YD* J ) ) GOTO 163 


0338 




166 


9TEMP=DTEMP*( 1.0-DR ELI J) ) 


0339. 




167 


CONTINUE — ~ 


0340 






DSUM=DSUM+DTEMP 


0341 




168 


CONTINUE 


0342 






WRITE! 6*961 DSUM 


0343 






GO TO 121 


0344 




163 


OTEMP=DTEMP*DREL(J ) 


0345 






GO TO 167 


0346 




165 


I F ( QN 1 ( YO t J ) ) GO TO 163 


0347 


r 




GO TO 166 - 




c 

r 




BLOCK NUMBER 27 


0348 




14 


DTEMP=DTEMP*DREL ( J ) 


0349 






GO TO 16 __ 


0350 




15 


DSUM=0 . 0 


0351 






USTP=URUN( R) 


0352 






DO 17— UO=URAN»USTP 


0353 






DTEMP=1 .0 


0354 






DO 16 J = 1 » E 


0355 






I F ( QPC ( UD » J ) ) C-0 TO 14 


0356 






OTEMP=OTEMP*( l.O-ORFLI J ) ) 


0357 




16 


CONTINUE 


0358 






DSUM=DSUM+DTEMP 


0359 




17 


CONTINUE 


0360 






WRITE <6, 98) R » ST ( P ) ? SP(R 


0361 






WRITE(6,95) OSUM 


0362 






WRITE(6»85) 


0363 






IF( SFL4.LE.0) GO TO 31 


0 364 






SELX=SEL4 


0365 




30 


I F ( SELX. GT . 15 ) GO TO 3 


0366 






DTEM1=DREL( SELX) 


0367 






OREL! SELX)=C.O 


0368 






0SUM=0.0 
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0369 




10 32 UD = URAN » USTP 


0370 




IF ( QPC ( U0* SELX ) ) GO TO 32 


0371- 




'lT(rHP=l * 0 


0372 




10 33 J = 1 * E 


0373 




IF (QPC(UD, J) ) CO TO 36 


0 3 74 




)TFMP=DTeMP* ( 1 ,0-DREL ( J ) ) 


03 75 


33 


CONTINUE 


0376 




OSUM=DSUM+QTEMF 


0377 


32 


CONT INUE 


037fi 




a'R I TE ( 6 » 45 ) SELX 


0379 




WR I TF ( 6 1 95 1 OS UN 


0380 




wRITE(6»85) 


0381 




1REL ( SELX ) =DTF N1 


0382 


31 


IF ( SEL5.LE.0 > GO TO 34 


0 383 




SFL X = SFL5 


0384 




SF L 5=0 


0385 




GO TO 30 


0 3 86 


36 


1TEMP=0TFMP*0RFL(J ) 


0387 




GO TO 33 


0388 


3 4 


IMSFL6.LE.0) GO TO 35 


0 3 89 




StLX=SEL6 


0390 




SFL6=0 


0391 




GO TO 30 


0392 - - 


33 


IF ( SEL7 . LE »0 ) GO TO 119 


0393 




IF( SFL7.GT. 15) GO TO 3 


0394 




0TEM1=0REL( SEL7) 


0395 




0RFL1 SHL7)=DCHG 


0396 




OSUM=0.0 


0397 




or 44 UO=URAN» USTP 


0398 




OTFMP=1.0 


0399 




00 43 J*1,E 


94-00 




I F < QPC < UD* J ) ) GO TO 42 


0401 




’)TFMP=DTFMP*< 1 . 0 - 0 R E L ( J ) ) 


04G2 




CONTINUE 


0403 




•)SUM=OSUM'HOTFMP 


0494 


4 4 


CONTINUE 


0 4C5 




kk I T F ( 6,46) SFL7,0CHG 


94C6 




.'RITE (6,95) DSIIN 


0487 




WR ITF (6,85) 


0408 




OR EL-( SFL7 )=0TEM1 


0409 




GO TO 119 


0410 


42 


OTEMP=DTEMP*ORFL(J ) 


0411 


C 


GO TO 43 


0412 


C 


t-ND 
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C GR = TOTAL NO OF PATHS AT START OF RUN 
C T IS USED AS A RUNNING INDEX 
C GNO(R) IS USED AS A TOTAL FOR THE RUN 
C IF SEL IS ZERO, CCMPUTER WILL HALT AFTER 
C THIS RUN; IF NONZERO, THEN AN ADDITIONAL 
C DATA CARD WILL BE READ AFTER THIS RUN 
C IF SFL 2 IS ZERO, NO Q-SUB-K FUNCTION 
C WILL BF FOUND; IF NONZERO, THEN Q-SU6-N 
C WILL BE DETERMINED 
C NOW COMPUTE Q-SUB-N 

C UNO(R) CONTAINS LENGHT OF Q-SUB-I,J(P) 

C R IS THE RUN NUMBER 

C ST IS THE STARTING TERMINAL 

C SP IS THE ENDING TERMINAL 

C W IS THE NO. OF INDEPENDENT CIRCUITS 

C K IS THE NO. OF POSS I BfcE- COMB INAT I ONS 

C OF THF GIVEN CIRCUITS 

C GR IS US FD AS A RUN-TOTAL-FR 

C ST 0° 99 MEANS THF NUMBER OF NODFS, EDGES, OR CIRCUITS 
C IS GREATER THAN THE LIMITS OF 15, 10, AND 10 
C STQO 999 MEANS THE LIMIT IF Q-SUB-K HAS BEEN REACHED 

C STOP 888 MEANS TOG MANY TERMS HAVE BEEN PRODUCED 

C AND ALL RESULTS ARE NOW IN DOUBT 
C SFL IS THE SELECTION CODE 

C IF SEL1 = 0, NO BFANCH RELIABILITIES ARE READ 
C IF SEL1 NOT ZERO, THE VALUE FOR Q-SUBI,J IS FOUND 
C IF SFL 2 = 2, THF VALUE FOR Q-SUB-N IS FOUND 
C MATRIX X BF CONTAINS THE INDEPENDENT CIRCUITS 
C MATRIX BCM IS XBF AS LOGICAL ELEMENTS 
C MATRIX XPTH CONTAINS THE GIVEN PATH SET 
C MATRIX BFC CONTAINS ALL POSSIBLE COMBIN- 
C ATI ON S OF BCM 

C MATRIX BP P CONTAINS THE RINGSUM OF THE 
C PATH AND BFC 

C MATRIX PFC IS THE PATH MATRIX, INDEXED 
C BY T OVER ALL RUNS 

C MATRIX XPFCI IS THE INTEGER VERSION OF PFC 
C MATRIX XY COUNTS THF NUMBER OF K-EOGED PATHS 
C MATRIX ENO CONTAINS BRANCH LOADING 
C MATRIX GNO HAS THF NUMBER OF PATHS 
C MATRIX 

C THIS ADDS UNIQUE PATH PRODUCT TERMS 
C TO THE COLLECTION NOW IN QPC 
C BUILD THE PATH PRODUCT TERMS 
C STORE FOR LATER TESTING OF UNIQUENESS 
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APPENDIX II 



The problems solved in examples (3.1-3) and (4.5) are combined 
and solved by the digital computer program. The output of the running 
of the program is included as Appendix II. 



THERE ARE 3 LINKS A_ND 3 NODES 

THE GIVEN CIRCUIT MATRIX IS 
t 1 I 



THE BRANCH RELIABILITIES ARE 

— u.9coocoo.9uoaaao-.Aaoc.co 



THE GIVEN 


PATH WAS 


















EOGF 1 


2 


3 4 5 6 


7 


B 


9 


10 


li 


12 


13 14 


15 


1 


0 


0 0 0 — 7) 




0 


0 


0 


0 


C 


0 0 


0 


MATRIX 


QPC 


FOLLOWS 


















EDGE l 


2 


3 4 5 6 


7 


8 


9 


To” 


11 


12 


13 14 


15 


— 1 
1 

1 

1 

0 


0 

1 

o 


0 

Q 

1 

1 

1 


— — 
















1 

1 



















— THER E ARE 2 PATHS 
THE PATH MATRIX IS 

EDGE l 2 3 4 5 6 7 8 9 10 11 12 13 14 15 

I o' *0 

C 1 1 

IN PAIRS, THE NUMBER OF K-EDGEO PATHS AND BRANCH LOADING 
EOGE l" 2 3 4 5 6 ~7 S' 9 10 ff 1 2 ~l 3 1~4 15 

111 
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2 



RUN NO. 1 FROM 1 TO 

Q-SUB-I , J = C. 98 09998 

THE VALUE OF 0- SUB'-TTTJ )~W ITHOUT BRANCH - F 
_ _ o=siiB-i , j = o . a aaaaaa_ — . — *. 



THE VALUE OF Q-SU8-(l,J) WITHOUT BRANCH 2 
0- SUB-1 , J” =~ 0 . 399999*9”’ 



THE VALUE OF Q-SUB-(I,J) WITHOUT BRANCH 3 
Q-SUB-I , J = 0.3999999 



THE VA4.UE_OF_ O-SUa^U-r J4— WITH BRANCH l 

HAVING RELIABILITY VALUEO . 500000 

Q-SUtt- I,J * O'. 9049999 
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THE 


GIVEN 


PATH 


WAS 


EDGE 


1 2 


3 A 


5 6 7 8 9 10 11 12 13 14 15 



t> i o~o a — cr — o — a ~o o o c o 'O cr 



MATRIX QPC 


FOLLOWS 
















EDGE 1 2 


3 


A 


5 6 7 


8 9 


10 


11 


12 


13 


14 


15 


• 0 1 


0 




















1 1 


C 




















0 1 


I 




















1 1 


I 




















1 0 


1 




















THERE ARE 


2 


PATHS 


- ' — 










’"'v” 


/ 


THE. PATH MATRIX 


IS . 


. . 








— 


— . — 


- — 


EDGE 1 2 


3 


4 


5 6 7 


8 9 


10 


1 I 


12 


13 


14 


15 


0 1 


0 




















1 0 


1 




















IN PAIRS, 


THE 


NUMBER OF K 


-EDGED 


PATHS 


AND 


BRANCH LOADING 


EDGE”! 2 
1 1 


•x 

-0- - 


4 


5 6 7 


8 9 


1C 


11 


12 


13 


14 


15 



1 l l 
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RUN NO 



1 TO 



3 



2 FROM 

Q-SUB- 1 t J = 0.93C9993 

THE - VALUE OF Q-SUB- TTTT)~WITHOUT BFANCH "1 
.0=£li&?.l.J = 0.89999.95 



THE VALUF OF Q-SUB-(ItJ) WITHOUT BPANCH 2 
Q-SUB- I * J = 0 . ft 0 9 9*9 9'9 



THE VALUE OF Q-SUB-(I,J) WITHOUT BRANCH 3 
Q-SUB- I » J = C. 3999999 



THE- -VALUE- OF Q-SUB-44~» J ) WITH oRANCH 2 

HAVING RELIABILITY VALUFC . 50000'* 

Q-SUB- I, J * 0 .9049999 
MATRIX QK FOLLOWS 

- -4 1- C- - 

1 0 1 
111 
0 1 1 

Q-S UB- K = 0.9719 998 
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THE GIVEN PATH WAS 

EDGE 1 2 3 A 5 6 7 ft 9 10 11 12 13 14 15 

0 0 1 0 C 0 0 0 c 0 0 C 0 0 0 

MATRIX QPC FOLLOWS 

FDGE 1 2 3 4 3 6 7 8 9 10 1 1 1 2 1 3 14 1 5 

- c- 0 1 

1 0 1 

Oil 

1 1 1 

1 1 0 



THERE ARF 2 PATHS 
THE PATH MATRIX IS 

EDGE 1 2 3 4 8 6 7 8 9 1? 11 12 13 14 15 

C 0 1 

1 1 G 

IN D A I R S » THF NUMBER OF K-EDGED PATHS AND BRANCH LOADING 
EDGE 1 "2 34 56 7 8 Q 10 11 12 13 H 15 

1 1 1 
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/ 



3 f £ i } M 



3 



RUN NO. 


3 F R UM 2 TU 3 


q-sub-i , j 


= C. 9809993 


THE VALUE 


OF O-SUB- ( I , J J WITHOUT BRANCH 


CKSLLB-I t J 


= C. 89999.99 


THE VALUE 


OF Q— SUB— ( I » J ) WITHOUT BRANCH 


O-SUB- I , J 


= 0 . 8 9° 999 9' 


THE VALUE 


OF Q-SUB-(I,J) WITHOUT BPANCH 


Q-SUB-I ♦ J 


- C. 8099999 


THE VALUE 


OF O-SUB- ( I i J I WITH BRANCH 



HAVING RELIABILITY VALUEO . 50000° 



Q-SUB-I » J 


= 0 . 9^49999 


MATRIX QK 
1 0 
0 1 
1 1 
I 1 


FOLLOWS 

I 

1 

1 

0 


Q-SUB-K = 


0.9719998 



83 



THE TOTAL 


NUMBER OF PATHS FOUND IN EACH 


RUN 


1 c. ""■ 1 a. 


t: 








IN PAIRSt 


THE 


Number of 


K-EDGED PATHS AND 


BRANCH LOADING 


RUN NU. 




I F POM 


1 TO 2 


- 


1 


1 


0 






1 


I 


i 












11 1 




RUN NU. 




2 FROM 


<*> 

! 

O 

h- 

9—4 




1 


1 


0 




- 


I 


1 


1 


— — — 




RUN NO. 




3 FROM 


2 TO 3 




I 


1 


C 






1 


1 


I 






TOTAL BRANCH LOADING 







3 


3 


3 






- - — — 




— 


. .. _ - _ 


- ■ 
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